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ABSTRACT 

We report the first abundance analysis of 17 M giant stars in the inner Galactic 
bulge at (l,b) = (0°, — 1°), based on R = 25,000 infrared spectroscopy (1.5 — 
1.8 /mi) using NIRSPEC at the Keck telescope. Based on their luminosities and 
radial velocities, we identify these stars with a stellar population older than 1 
Gyr. We find the iron abundance ([Fe/H]) = —0.22 ± 0.03, with a la dispersion 
of 0.14±0.024. We also find the bulge stars have enhanced [alpha/Fe] abundance 
ratio at the level of +0.3 dex relative to the Solar stars, and low < 12 C/ 13 C >~ 
6.5±0.3. The derived iron abundance and composition for this inner bulge sample 
is indistinguishable from that of a sample of M giants our team has previously 
studied in Baade's Window (I, b) = (0.9, —4). We find no evidence of any major 
iron abundance or abundance ratio gradient between this inner field and Baade's 
window. 

Subject headings: Galaxy: bulge — Galaxy: abundances — stars: abundances 
— stars: late-type — techniques: spectroscopic — infrared: stars 



1 Data presented herein were obtained at the W.M.Keck Observatory, which is operated as a scientific 
partnership among the California Institute of Technology, the University of California, and the National 
Aeronautics and Space Administration. The Observatory was made possible by the generous financial 
support of the W.M. Keck Foundation. 
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Introduction 



High resolution infrared (IR) spectroscopy now enables the possibility of undertaking 
abundance analysis for cool late-type giants. The technique is especially powerful in ex- 
tending abundance studies to regions of high extinction, such as the inner disk and bulge. 
We h ave reported the first detailed abundances for Galactic bulge M giants (IRich fc Origlia 
20051 ) in Baade's Window, a field with low enough extinction that optical high resolution 



spectroscopy is also feas ible (IMc William fc Rich! Il994j ; iFulbright. McWilliam fc Richl 12006 



20071 ; IZoccali et al.l I2006T ) . Our IR spectroscopy has been in basically good agreement with 
these optical studies, finding roughly Solar mean abundance and enhanced alpha elements, 
consistent with a predominant chemical enri c hment by type II SNe on a relatively short 



timescale (IMatteucci. Romano &: Molarolll999l ; iBallero et al. 



20071 ). 



Here we exploit these techniques to push into the high extinction regions of the inner 
Galactic bulge, where optical spectroscopy is impossible. Our goal is to constrain whether 
an abundance or abundance r atio gradient is p resent in the bulge population. The Galactic 
Center contains old stars (e.g. iFiger et al.ll2004i ) and th e bulge in Baade's Window and other 
fields has been demonstra ted to be globular cluster age (jOrtolani et al.lll995l ; iKuijken fc Rich 
20021 ; IZoccali et al.l 120031 ). Our aim has been to select the old red giant stars that lies near 
the red giant branch (RGB) tip a factor of four closer to the nucleus than the well studied 
Baade's window field, at (l,b) = (0°, — 1°), i.e. only 140 pc distant from the nucleus 
(adopting R = 8 kpc), with the goal of measuring possible abundance gradients in the inner 
Galactic bulge. 



At present, there is no secure measurement of the bulge abundance gradient. iFrogel et al 



( 119991 ) used broadband JK photometry to probe fields along the major and minor axis from 
—4° < l,b < 0°, reaching within 0.2° of the Galactic Center. This study, which found n o 
indication of a significant gradient, provided the input sample for iRarmrez et al.l (120001 ). 
who also found no gradien t using low resolution measurements of Ca, Na, and CO features 
in the 2/zm window. The iRamfrez et al.l (120001 ) abundances were derived from a globular 
cluster-based calibration and consequently their abundances for [Fe/H]>0 are an extrapola- 



tion; h owever, that was unlikely to have altered their findings. More recently, IZoccali et al. 
( 120031 ) find a photometric abundance distri bution at (I, b) = (0.3 °, —6.2°) ve r y sim ilar to 



that of IFulbright. McWilliam fc Richl (120061 ) in Baade's window. IViera et al.l (120061 ) finds 
the same photometric abundance distribution in a proper motion selected sample of K gi- 
ants in the Plaut field, at b = —8°. These studies suggest that the bulge does not have a 
measurable abundance gradient over its inner kpc, but they all rely on photometric abun- 
dance constraints (RGB color) rather than on spectroscopy. 
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Fig. 1. — NIRSPEC H-band spectrum showing the region near 1.555 fim of three (#9, #10 
and #11) of the coolest and more metal rich stars in our sample. A few major atomic lines 
and molecular bands of interest are flagged. 
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2. Observations and abundance analysis 



A sample of M giants near the RGB tip in the inner Bulge field at (/, b) = (0°, —1°) has 
been selected from the 2MASS K, J-K co lor magnitude diagram , corrected for reddening by 
using the interstellar extinction maps of ISchultheis et al.l (119991 ) . Spectra for 17 stars with 
bolo metric magnitudes Mboi > —4 were ta ken during 3 observ ing runs on April 2005, May 
2005 and May 2006, by using NIRSPEC jMcLean et al.l Il998h at Keck II. A slit width of 
0'/43 giving an overall spectral resolution R=25,000, and the standard NIRSPEC-5 setting, 
which covers most of the 1.5-1.8 yum H-band have been selected. 

The raw stellar spectra have been reduced using the REDSPEC IDL-based package 
written at the UCLA IR Laboratory. Each order has been sky subtracted by using nodding 
pairs and flat-field corrected. Wavelength calibration has been performed using arc lamps 
and a second order polynomial solution, while telluric features have been removed by using 
a O-star featureless spectrum. The signal to noise ratio of the final spectra is >30. Fig. [T] 
shows an example of the observed spectra for 3 among the coolest and most metal rich giants 
in our sample. 

A grid of suitable synthetic spectra of giant stars has been computed by varying the 
photospheric parameters and the element abun dances, using an updated version of the code 
described in lOriglia. Moorwood fe Olival (119931 ). By combining full spectral synthesis anal- 
ysis with equivalent widths measurements of selected lines, we derive abundances for Fe C, 
O and ot her alpha-elements (Mg, Si, C a and Ti). The l i nes a nd analysis method are de- 
scribed in lOriglia. Rich fe Castrol (120021 ); lOriglia fe Richl (120041 ) and subjec ted to rigorous 
tests in our previous studi es of Galactic bulge field and cluster giants (see lOriglia fe Rich 



2004 



Rich fe Orig 



Grevesse fe Sauval 



laM/iUUo 
Jl998h . 



2005, and references therein). Reference solar abundances are from 



In the first iteration, we estimate stellar temper ature from the (J — K ) n colors (see 
Tabled]) by using an average red ddening of E(B-V)=2.9 (Schultheis et aU1999l ) and the color- 
temperature transformation of iMontegriffo et al.l (119981 ) specifically calibrated on globular 
cluster giants. Gravity has been estimated from theoretical evolut ionary tracks, acco rding 
to the location of the stars on the Red Giant Branch (RGB) (see lOriglia et al.lll997l . and 
references therein for a more detailed di scussion). An avera ge value £=2.0 km/s has been 
adopted for the microturbulence (see also lOriglia et al.lll997l ). More stringent constraints on 
the stellar parameters are obtained by the simultaneous spectral fitting of the several CO and 
OH molecular bands, whic h are very sensitive to t e mper ature, gravity and microturbulence 
variations (see Figs. 6,7 of lOriglia. Rich fe Castrol (120021 )). 



CO and OH in particular, are extremely sensitive to T e ff in the range 3500 to 4500 
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K. Indeed, temperature sets the fraction of molecular versus atomic carbon and oxygen. 
At temperatures >4500 K molecules barely survive; most of the carbon and oxygen are in 
atomic form and the CO and OH spectral features become very weak. On the contrary, at 
temperatures <3500 K most of the carbon and oxygen are in molecular form, drastically 
reducing the depe ndence of the CO a nd OH band strengths and equivalent widths on the 



temperature itself (jOriglia et al.lll997l ). 



The final values of our best-fit abundances together with random errors are listed in 
Table [U 

As a further check on the statistical significance of our best-fit solution, we also compute 
synthetic spectra with AT c j = ±200 K, Alog g = ±0.5 dex and A£ = =(=0.5 km s _1 , and with 
corresponding simultaneous variations of the C and O abundances (on average, ±0.2 dex) to 
reproduce the depth of the molecular features. As a figure of merit of the statistical test, we 
adopt the difference between the model and the observed spectrum (hereafter 5) . In order to 
quantify systematic discrepancies, this parameter is more powerful than the cla ssical y 2 test. 



which is instead equally se nsitive to random and systematic scatters (see also lOriglia et al. 



20031 ; lOriglia fc Richl 12004 ) . Our best fit solutions always show >90% probability to be 
representative of the observed spectra, while those with ±0.1 dex are significant at >1 a level. 
Spectral fitting solutions with abundance variations of ±0.2 dex, due to possible systematic 
uncertainties of ±200 K in temperature, ±0.5 dex in gravity or =(=0.5 km/s in microturbulence 
have <30% probability of being statistically significant. Hence, as a conservative estimate 
of the systematic error in the derived best-fit abundances, due to the residual uncertainty in 
the adopted stellar parameters, one can assume a value of < ±0.1 dex. However, it must be 
noted that since the stellar features under consideration show a similar trend with variations 
in the stellar parameters, although with different sensitivities, relative abundances are less 
dependent on the adopted stellar parameters (i.e. on the systematic errors) and their values 
are well constrained down to ~ ±0.1 dex (see also Tabled]). 



3. Results and Discussion 



For the observed M giants we measure average < [Fe/H] >= — 0.22±0.03 and dispersion 
la = 0.140 ± 0.024. Th is [Fe/H| abundance distribution is very similar to the one measured 
in the Baade's window (IRich fc Orighall2005l ). with almost identical average value within the 
error and only marginally greater dispersion. 



We also confirm the lack of both meta l poor and super metal rich M giants. This was 
originally noted in the lRich fc Origlial (120051 ) Baade's Window M giant distribution, but now 
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Table 1. Stellar parameters and abundances for our sample of giant stars in the inner 

Bulge field at (l,b) = (0,-1). 



it 


RA 




Dec 


(J - K1S 


T a 
A en 






fFe/Hl 


fO /Fel 


fSi/Fcl 


fMe /Fel 


[Ca/Fe] 


fTi/Fel 

L / J 


fC /Fel 


i 


17:49:53.6 


-29 


:31:22.4 


0.85 


4000 


1.0 


-327 


+0.02 


+0.28 


+0.28 


+0.30 


+0.28 


+0.38 


-0.32 


















± 0.11 


± 0.16 


± 0.20 


± 0.16 


± 0.19 


± 0.20 


± 0.17 


2 


17:49:43.2 


-29: 


:18:50.9 


1.25 


3400 


0.5 


-66 


-0.40 


+0.41 


+0.30 


+0.32 


+0.30 


+0.23 


-0.40 


















± 0.08 


± 0.11 


± 0.21 


± 0.09 


± 0.12 


± 0.14 


± 0.11 


3 


17:49:53.6 


-29: 


:30:41.2 


1.04 


3600 


0.5 


-76 


-0.23 


+0.22 


+0.33 


+0.32 


+0.33 


+0.33 


-0.22 


















± 0.08 


± 0.11 


± 0.17 


± 0.08 


± 0.12 


± 0.14 


± 0.10 


4 


17:49:46.5 


-29. 


19:17.8 


1.37 


3200 


0.5 


-7 


-0.11 


+0.32 


+0.21 


+0.30 


+0.31 


+0.22 


-0.29 


















± 0.07 


± 0.11 


± 0.16 


± 0.08 


± 0.12 


± 0.13 


± 0.10 


5 


17:49:51.2 


-29: 


:31:00.4 


1.21 


3400 


0.5 


-192 


-0.20 


+0.22 


+0.25 


+0.24 


+0.30 


+0.18 


-0.25 


















± 0.08 


± 0.12 


± 0.17 


± 0.09 


± 0.12 


± 0.15 


± 0.11 


6 


17:49:36.7 


-29: 


:18:53.8 


1.40 


3200 


0.5 


-32 


-0.40 


+0.27 


+0.24 


+0.27 


+0.30 


+0.25 


-0.30 


















± 0.08 


± 0.11 


± 0.14 


± 0.09 


± 0.12 


± 0.14 


± 0.11 


7 


17:49:40.8 


-29: 


:18:17.2 


1.09 


3600 


0.5 


+45 


-0.24 


+0.30 


+0.34 


+0.34 


+0.34 


+0.34 


-0.26 


















± 0.03 


± 0.05 


± 0.13 


± 0.04 


± 0.05 


± 0.07 


± 0.08 


8 


17:49:53.3 


-29: 


30:28.02 


1.10 


3600 


0.5 


-1 


-0.17 


+0.26 


+0.27 


+0.29 


+0.27 


+0.35 


-0.33 


















± 0.03 


± 0.05 


± 0.13 


± 0.04 


± 0.05 


± 0.07 


± 0.08 


9 


17:49:47.0 


-29: 


:18:49.2 


1.47 


3200 


0.5 


-15 


-0.10 


+0.27 


+0.22 


+0.29 


+0.30 


+0.23 


-0.30 


















± 0.04 


± 0.06 


± 0.13 


± 0.04 


± 0.06 


± 0.07 


± 0.08 


10 


17:49:47.6 


-29. 


:18:59.8 


1.42 


3200 


0.5 


+ 121 


-0.17 


+0.19 


+0.22 


+0.26 


+0.27 


+0.27 


-0.23 


















± 0.04 


± 0.06 


± 0.13 


± 0.04 


± 0.06 


± 0.07 


± 0.08 


11 


17:49:35.5 


-29: 


18:46.1 


1.49 


3200 


0.5 


-168 


-0.07 


+0.15 


+0.17 


+0.24 


+0.27 


+0.22 


-0.33 


















± 0.04 


± 0.06 


± 0.13 


± 0.04 


± 0.06 


± 0.07 


± 0.08 


12 


17:49:41.8 


-29 


:18:25.4 


1.36 


3200 


0.5 


-94 


-0.10 


+0.16 


+0.20 


+0.29 


+0.30 


+0.30 


-0.25 


















± 0.09 


± 0.14 


± 0.17 


± 0.10 


± 0.13 


± 0.15 


± 0.11 


13 


17:49:39.5 


-29 


20:12.3 


1.00 


3800 


0.5 


+ 165 


-0.34 


+0.43 


+0.44 


+0.41 


+0.34 


+0.40 


-0.26 


















± 0.09 


± 0.10 


± 0.21 


± 0.10 


± 0.13 


± 0.14 


± 0.12 


14 


17:49:49.0 


-29: 


19:50.8 


1.07 


3400 


0.5 


-181 


-0.11 


+0.27 


+0.26 


+0.27 


+0.31 


+0.19 


-0.29 


















± 0.07 


± 0.11 


± 0.19 


± 0.08 


± 0.12 


± 0.15 


± 0.10 


15 


17:49:39.4 


-29: 


:20:21.3 


1.33 


3400 


0.5 


-34 


-0.42 


0.37 


0.32 


0.32 


0.32 


0.26 


-0.38 


















± 0.08 


± 0.11 


± 0.15 


± 0.09 


± 0.12 


± 0.14 


± 0.11 


16 


17:49:40.7 


-29: 


:20:13.3 


1.35 


3600 


0.5 


-90 


-0.46 


0.37 


0.36 


0.34 


0.36 


0.36 


-0.24 


















± 0.08 


± 0.10 


± 0.15 


± 0.09 


± 0.12 


± 0.14 


± 0.11 


17 


17:49:35.4 


-29: 


:17:45.0 


1.25 


3400 


0.5 


-206 


-0.25 


+0.22 


+0.19 


+0.23 


+0.25 


+0.25 


-0.35 


















± 0.08 


± 0.11 


± 0.17 


± 0.09 


± 0.12 


± 0.14 


± 0.10 



12, 



a The (J-K) colors are from 2MASS and have been corrected for reddening using E(B-V)=2.9. 
b Hclioccntric radial velocity in km s — 1 . 
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Fig. 2. — Histogram of the metallicity distribution in [Fe/H], for the observed giants 
in the inner bulge field at (7, 6) = (0°, — 1°) (shad ed histogram) and in Baade's window 
(Z,6) = (0.9°, -4°) (solid line. iRich fc Oririial (booi l). for comparison. Notice that there is 
no significant difference between the two fields. 
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Fig. 3. — a/Fe as a function of [Fe/H] for the observed giants in the inner bulge field at 
(I, b) = (0°, —1°) (filled circles) and in Baade's window (open circles. iRich fc Origfial (120051 )). 
for comparison. Typical errors are plotted in the bottom-right corner of each panel. There 
is no evidence for a gradient in [a/Fe] in the bulge. 
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persists in a sample that, in total, is nearly three times as large. We know that the metal 
poor end will be poorly represented, as few stars with [Fe/H]< —1 will evolve to be cool 
enough to develop TiO bands and become M giants at the RGB tip. However, the lack 
of super metal rich M giants is a conc e rn. In deed, although they represent only a fraction 



(<20%, iFulbright. Mc William fc Rich! (120061 )) of the bulge population, they are definitely 



present in all bulge K giant studies. We speculate that their low intrinsic numbers and rapid 
evolution may make their detection statistic ally unlikely. Moreover, the possible enhanced 



mass loss at super solar metallicities (see e.g. ICastellani fc Castellanilll993l ) may also prevent 



the most metal rich stars from reaching the RGB tip. It is important to explore these issues 
with a larger sample, but the total of 31 M giants observed in these two fields is now large 
enough that the lack of stars with [Fe/H]>0.1 is a significant concern. 

Our survey shows a rather homogeneous a-enhancement by «+0.3±0.1 dex up to solar 
metallicity, without significant differences among the various a-elements (see Fig. [3]), and 
similar to measurements in the Baade's window. Such an enhancement of the [a/Fe] abun- 
dance ratio suggests that at the epoch of the bulge formation, mainly type II SNe were 
contributing to the enrichment of the interstellar medium. It also indicates that the overall 
formation process should have ended before the bulk of type la SNe have exploded, i.e. well 



withi n a few Gyr, the precise timescale depending on the SN modeling (see e.g. iGreggio 



20071 ). 



We have also derived carbon abundances for all of the observed stars from analysis of the 
CO bandheads. We find some degree of [C/Fe] depletion (< [C/Fe] >= -0.29 ± 0.02) with 
respect to solar values and very low 12 C/ 13 C (< 12 C/ 13 C >= 6.5 ± 0.3) , as also measured 
in Baade's window M giants, confirming the presence of extra-mixing processes during the 
RGB phase of evolution in metal-rich environments. 

From our detailed analysis of chemical abundances and abundance patterns we can thus 
exclude any major abundance gradient in the inner bulge. 

Finally, we have also measured radial velocities and velocity dispersion for the observed 
stars, as for the Baade's window giants. The overall velocity dispersion of the global sample 
of 31 M giants turns out to be cr Vr ~ 116 ± 15 km/s, fully consistent with bulge kinematics. 



Theory concerning the formation of spheroids has evolved, beginning with lEggen et al. 



(119621 ). who proposed violent relaxation and dissipative collapse. In the CDM scenario, 



spheroids would form from the mergers of disks. Further, large galaxy sur veys now sugges t 



that the spheroid-dominated " red sequence" has grown in mass since z ~ 1 (jBell et al.ll2004j ) , 
implying that some spheroids have come into existence relatively later than the age of the 
Galactic bulge population. A widely supported scenario for bulge formation was initially 
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proposed by iRaha et al.l (ll99ll ). In this scenario, a massive disk develops bending modes 
that ultimately thicken the disk vertically, producing a p eanut-shaped bar. These idea s have 
been developed by many other theoretical studies e.g. iPfenniger &: Norman! (Il990l ). The 
bar thickening hypothesis predicts that no vertical abundance gr a dient should be present. 
Carefully weighing the existing evidence, iKormendy fc Kennicuttl (120041 ) conclude that the 
Galactic bulge has likely been formed via secular processes, but they are perplexed by the 
strong evidenc e for its age bein g similar to that of the halo. The recent survey of bulge 
kinematics by iRich et al.l (120071 ) supports the idea that the bulge is kinematically distinct 
from the i nner disk a nd that while the deprojection of the bulge demands some elongated 
structure (lZhadll996l ). the new velocity field also requires the presence of retrograde orbits 
by construction. While the lack of an abundance gradient is consistent with secular evolution 
models, the great age of the bulge/bar remains as a problem for this picture. 
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